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Abstract 

We introduce an exactly-solvable family of one-dimensional driven- 
diffusive systems defined on a discrete lattice. We find the quadratic 
algebra of this family which has an infinite-dimensional representation. 
We discuss the phase diagram of the system in a couple of special cases. 

One-dimensional driven-diffusive systems are systems of classical particles 
with hard-core interactions moving in a preferred direction which can be used to 
model many systems such as ribosomes moving along a m-RNA, ions diffusing in 
a narrow channel under the influence of an electric field or even cars proceeding 
on a long road [TJ [SJ [3] . These systems are usually defined on an open lattice 
coupled with two reservoirs at both ends or on a lattice with periodic boundary 
conditions. In long time limit the system settles into an non-equilibrium steady 
state characterized by some bulk density and the corresponding particle current. 
The out-of-equilibrium steady state properties of these systems are not only 
affected by the boundaries but also to the localized inhomogeneities which play 
a crucial role not comparable to that in classical equilibrium systems. Driven- 
diffusive systems have been studied extensively in the past decade because of 
their unique non-equilibrium properties. 

There are different approaches to study the steady state of these systems. 
The Matrix Product Formalism (MPF) is one of the most powerful techniques 
in this field [H [H [5] . According to this formalism one assigns an operator to ev- 
ery state of a lattice site of the system. Now every configuration of the system 
is associated with a product of such operators. The steady state probability 
of such configuration is then given by a trace (for the systems defined on a 
ring geometry) or a matrix element (for the systems defined on an open lat- 
tice) of such products. For instance for a three-states system on a lattice with 
periodic boundary condition we define three operators, let us say A, B and 
E, associated with three different states of each lattice site. The steady state 
weight of a given configuration similar to ABEEABEAA is then proportional 
to Tr[ABEEABEAA]. Requiring that the probability distribution function of 
the system defined above is stationary provides us with an algebra of operators. 
For the systems with nearest neighbors interactions it turns out that the algebra 
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is quadratic. In order to calculate the steady state weights one can either work 
with the commutation relation of the operators or find a representation for the 
associated quadratic algebra of the system (for a recent review on the MPF see 

ED)- 

So far only a couple of one-dimensional driven-diffusive models have been 
solved exactly using the MPF 0]. Therefore it would be interesting to investi- 
gate whether or not there are other models which can be solved exactly using 
this technique. In this paper we introduce a large family of one-dimensional 
driven-diffusive models which under some constraints on the reaction rates can 
be solved exactly i.e. the steady state weight of any configuration can be cal- 
culated rigourously using the MPF. From there the mean values of physical 
quantities such as the currents of particles and also their concentrations can 
be obtained exactly. This three-states family belongs to the systems with non- 
conserving dynamics and nearest-neighbors interactions. These states belong 
to two different types of particles and holes on each lattice site. It is not clear 
that the steady state weights of such a system can generally be written as ma- 
trix product states; however, we will present a constraint under which it will 
possible to write these weights using the MPF. In this regard, we will relate the 
quadratic algebra of our system to the quadratic algebra of the systems with 
known representations. A couple of members of this family of systems have al- 
ready been studied in literature. At the end of this paper we will briefly discuss 
some other special cases which have not been studied before. 

During last decade several quadratic algebras have been introduced and used 
to study the critical behaviors of one-dimensional driven-diffusive systems. It 
has also been tried to classify certain quadratic algebras and find their represen- 
tations [7, 8, 9j. One of the most well-known quadratic algebras belongs to the 
Partially Asymmetric Simple Exclusion Process (PASEP) in which the classical 
particles hop to the left and to the right with the rates p and q on an open 
one-dimensional lattice of length L. The jumps of particles are only successful 
provided that the target sites are empty. The particles are injected into the 
system from the left boundary with the rate a provided that the first site of the 
lattice is empty. The particles are also extracted from the last site of the lattice 
provided that it is already occupied. The particles are subjected to the hard- 
core interactions so that two particles cannot occupy a single site of the lattice 
simultaneously. It is known that the steady state weights of the PASEP can 
be written in terms of products of infinite-dimensional square matrices which 
satisfy a quadratic algebra of the following form [3 [10l [11] 

pAB - qBA = A + B 

aA\V) = \V) (1) 
(3{W\B = (W\ 

in which the operators A and B are associated with two different states of each 
lattice site i.e. a particle and an empty site respectively. It is known that this 
algebra has an infinite-dimensional matrix representation for p > q given by the 
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following matrices |5j 
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and vectors 
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(W\ = ( 1 ••• ) 



(3) 



in which we have defined a = 1 



b = 1 



and ci- 



(l-(f) 2 )(l 



/3 ^p' 

1 ). The matrix representation of |T|) given by ([2]) and {3} is quit well- 
defined in a sense that product of any number of these matrices sandwiched 
between (W\ and \V) is finite. In what follows we will look for those systems 
which their steady states in terms of the MPF are given by ([T]) and its repre- 
sentation i.e. and 

Let us define a new operator E — u>\V)(W\ in which uj is a real number. By 
multiplying (W|u; from the right in the second row and also oj\V) from the left 
in the third row of ([T]) one finds 



pAB - qBA = A + B 
aAE = E 
(3EB = E 



(4) 



which was first introduced in |12j . Obviously this algebra has the same repre- 
sentation given by (|2|) and ((3]) and one should only define E as 

/ lu ■•• \ 









E = uj\V)(W\ 



\ 



(5) 



/ 



It can easily be verified that the operator E has the property E 2 



lE and 



also can be added to {4|) without changing its matrix representation to make 
a larger and well-defined quadratic algebra. Therefore the following quadratic 
algebra has the same matrix representation given by ^ and {5} 



pAB - qBA = A 
aAE = E 
(3EB = E 
E 2 = uE. 



B 



(6) 
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The question is now whether or not one can find a three-states model, defined 
on a lattice with a ring geometry, with a steady state described by By 
applying the standard MPF we have found that the steady state probability 
distribution function of the following three-states system 

A + B -> B + A with rate x 42 

B + A -> A + B with rate x 2 4 

A + -> + A with rate x 73 

A + % -> + B with rate ir 83 

A + -> + with rate £ 93 , . 

+ 0^^ + with rate a; 39 ^ ^ 

+ B -> A + with rate £ 38 

+ fl -> fl + with rate £ 68 

+ B -> + with rate £ 98 

+ -> + B with rate a; 89 

in which A, B and _E can be associated with two particles of different types and 
an empty site respectively, can be described by ^ provided that we define 

P = X42 

q = a; 24 

a = x 73 (8) 

[3 = x 68 

LU= — i (i£2i + 22a) 

and require the parameters to satisfy the following constraint 

,£83 X 38 + X 9S X 3S X 93 +X 83 

x 3 g{ ) = x S9 { ). (9) 

X73 x 68 x 6 8 x 73 

This constraint together with the definitions ([5]) guarantee that our model de- 
fined by ||7J| has a well-defined algebra given by © and that its steady state 
weights can be expressed in terms of a matrix product states. 

As we mentioned, a couple of special cases of this family of models have 
already been studied in literature. For instance in [13] the authors have studied 
an exactly solvable model with the constraints a; 83 = x 3 s — 0, X42 = ^24 = 
xj 3 = xqs — x 3 g — x$g = 1 and xg 3 — xqs = w. In this case the particle number 
for both A and B particles changes because of creation and annihilation of them. 
They have found that the system undergoes a continuous phase transition from 
a fluid phase to a maximal current phase for q < 1 by varying uj . For q > 1 the 
system is always phase separated and there is no phase transition. In another 
paper the authors have considered an special case in which X24 — x^ 3 = x 3 g = 
x 9 8 = .r 89 = 0, X42 = xqs = £73 = x 3 q = 1 and xg 3 = to [H]. In this case 
the number of B particles is conserved while the number of A particles can 
change due to the creation and annihilation. It has been shown that, in this 
case where a finite density of B particles exists on the ring, the model can be 
solved exactly and undergoes a second-order phase transition by varying u) from 
a phase in which the density of the empty sites is zero to another phase where the 
density of the empty sites is nonzero. In |15j the authors have studied the same 
process with more general reaction rates as X24 — £ 83 = £ 38 = x 98 = a; 89 = 0, 
X42 = 7^68 = a X73 = X39 ~ 1 an< ^ X93 = w ' nowever j m this case they assume 
that there is only a single B particle in the system. They have found that 
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the phase transition is now discontinuous and that shocks might appear in the 
system at the transition point. 

In what follows we consider the most general case where most of the pa- 
rameters in ([7]) are nonzero which has not been studied in literature yet. We 
first study the phase diagram of the model for X24 — and X42 = 10- The 
partition function of the system defined as sum of the weights of all accessible 
configurations with at least one empty site, is given by 

Z L (a,f3,u,0 = Tr[(£A + B + E) L ]-Tr[(£A + B) L ]. (10) 

The fugacity of A particles £ is an auxiliary parameter and the reason we have 
defined it will be clear shortly. It turns out that the generating function of 
this partition function can be calculated exactly. After some straightforward 
calculations one finds 

G(a,/?,u;,£,A)= VA^^-^- (11) 
l=i 1 _ uU 

in which we have defined 

U(a,P,t,\) = J2 xL+1 (Wm + B) L \V). 

It turns out that 



(12) 



where 



1 



/±0r) = -(1 - 2x ± A(l - - V(l + A(l-£)) 2 -4A). 



.7: 



The phase diagram of the system can now be obtained by studying the sin- 
gularities of the generating function (JTTJ) for £ = 1. One can easily see that 
in this case the generating function has two different kinds of singularities: a 
square root singularity A* = \ and two simple pole singularities which come 
from denominator of (jTTJ) by solving the equation 1 — LjU(a,(3,£ — 1, A*) = 0. 
However, analyzing the absolute values of singularities shows that the system 
can only have two phases: a maximal current phase which is specified by the 
square root singularity (for a + f3 > 1 and (2 — ~) (2 — 4) > lo) and a fluid phase 
which is specified by the simple pole singularity (for a + (3 < 1 or a + /3 > 1 and 
(2 — i)(2 — i) < uj). Keeping £ = 1 the total density of the empty sites pe can 
easily be calculated using 

p E = lim ^^-lnZ L . (13) 

L^oo L OU! 

In the thermodynamic limit we have pe ~ — cj 31 ^ . It turns out that the 
density of the empty sites is zero in the maximal current phase while it is 
nonzero in the fluid phase. The total density of A and B particles in each phase 
can also be calculated exactly. Since the density of the empty sites is known 



x By resealing the time one can always take one of the parameters equal to one. 
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only one of these densities is independent. The auxiliary fugacity £ can now 
help us find the density of A particles as 

The density of the B particles is then ps = 1 — Pe ~ PA- We have found that 
in the maximal current phase pa = Pb = \- In the fluid phase both pa and ps 
are complicated functions of a, (3 and u> and will not be presented here. The 
particle currents for both species can also be calculated exactly in this case. 
We have found that the current of A particles J a is always equal to that of 
B particles Jb- Our calculations also show that the particle current is given 
by J a = Jb = Z 'z L 1 which is equal to A* in the thermodynamic limit. In the 
maximal current phase we simply find them to be equal to \. This is actually 
in contrast with the case studied in [14] where the number of B particles in one 
species is conserved. It has been shown that in this case the currents can be 
different. For X24 = q, 2:42 = 1, a = 1 and (3=1 the results are exactly those 
obtained in [T3] . On a lattice with periodic boundary conditions we have found 
that the phase diagram of the model does not change even for arbitrary a and 
/?; therefore, we will not discuss this case here. 

One should note that ([6]) can also explain the steady state of a system with 
open boundaries and two species of particles. The particles of type A (B) are 
injected from the left (right) boundary with rate — (— ) and extracted from the 
right (left) boundary with rate a (/?). All of the processes in (J7J) might also 
take place on the lattice. The partition function of the model for 224 = and 
£42 = 1 can also be calculated exactly and is given by 

Z L (a, M = (W\(A + B + Ef\V) = ± ^I^^IP 1 (15) 
in which 



and 

Ai = — !— (a - (3 + af3u - J (a -(3 + af3uj) 2 + 4af3 2 co(l - a)), (17) 
2a[3uj 

A 2 = — —fa -(3- af3uj - J (a -(3 + af3uj) 2 + 4a/3 2 w( 1 - <*))• (18) 
2apuj 

The phase diagram structure of the model can now be obtained by studying 
the thermodynamic behavior of the partition function. Equivalently one can 
study the zeros of this partition function as a function of a, (3 or to. It turns out 
that the system has again two different phases. The phase transition occurs for 
a + (3>latuj — (2 — i)(2 — ^) similar to the ring geometry case. Both phases 
are symmetric that is the currents of particles of different types are always equal. 

In this paper we have investigated a general quadratic algebra §E§ associ- 
ated with a family of exactly solvable three-states reaction-diffusion systems 
with non-conserving dynamics defined on a one-dimensional lattice with ring 
geometry. This algebra is in fact the quadratic algebra associated with the 
PASEP given by |T]) besides the relation E 2 = u>E which does not change the 
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1 + auj\i 



1 + (3uj\ 2 



(16) 



representation of the algebra but it generates a new algebra which allows us 
to study a new family of three-states processes using the MPF. This family of 
three-states processes are defined by ten nonzero reaction rates given by ([7]) 
which should satisfy a constraint given by ©. Under this constraint the steady 
state of the system can be written as a matrix product form. We have considered 
the most general model of this type on a lattice with periodic boundary condi- 
tions and studied its phase diagram. As we have also mentioned, the generalized 
algebra ^ can explain the steady state of a three-states system with the reac- 
tions defined by ([7]) but this time under open boundary conditions where the 
particles are allowed to enter and leave the lattice with some specific injection 
and extraction rates. The phase diagram of the model has also been studied 
and the partition function calculated exactly. Our approach can be generalized 
and applied to other models similar to the model studied in [TB] (which is a 
p-species model defined on a lattice with periodic boundary conditions) to find 
a, p + 1-species exactly solvable model. The results will be published elsewhere. 
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